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Abstract. We review the formalism of quark recombination applied to the 
hadronization of a quark gluon plasma. Evidence in favor of the quark 
recombination model is outlined. Recent work on parton correlations, leading to 
detectable correlations between hadrons, is discussed. Hot spots from completely 
quenched jets are a likely source of such correlations which appear to be jet-likc. 
It will be discussed how such a picture compares with measurement of associated 
hadron yields at RHIC. 
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In the first part of this manuscript we revisit some basic arguments leading to 
the formulation of recombination. Applications to RHIC are highlighted, in particular 
the discovery of the quark counting rule for elliptic flow. Toward the end of the first 
part we discuss some topics fundamental to recombination that have been disputed 
recently. We refer to and the original papers cited therein for a more thorough 
review. 

In the second part we introduce a recent extension of the recombination formalism 
that permits correlations between partons [2]. These naturally lead to correlations 
between hadrons. We argue how jet-like correlations at intermediate Pt can be 
understood as signals of quenched jets. 



1. Recombination of Quarks: a Review 

1.1. Introduction 

In high energy collisions of particles or nuclei, partons are scattered out of the 
bound states and new partons are created. The confinement property of quantum 
chromodynamics (QCD) does not permit the existence of these colored partons as 
freely propagating states. Instead they have to be grouped into color singlet bound 
states that eventually become hadrons. This hadronization is a non-perturbative 
process in QCD and its dynamics is not fully understood. 

It has been proven that for a perturbative process at asymptotically large 
momentum transfer, the non-perturbative hadronization process can be separated 
from the hard scattering amplitude in a well-defined way. Hadronization in this limit 



Hadron Correlations from Recombination and Fragmentation 



2 



is described by a set of fragmentation functions D a i^ that parametrize, in a universal 
way, the probability that a hadron A with momentum P is created from a parton a 
with momentum p in the vacuum [3j . Fragmentation functions have been measured in 
e + + e~ collisions and work well for hadron production at transverse momenta Pt > 1 
GeV/c in p + p at RHIC energies. 

However, the fragmentation process is not sufficient to explain hadron production 
in Au+Au collisions at RHIC at transverse momenta of several GeV/c. Two puzzling 
observations lead to this conclusion. First, baryons are much more abundant than 
predicted by fragmentation. PHENIX has observed a ratio of proton/pion as 1 between 
2 and 4 GeV/c 0|. The value predicted by leading twist perturbative QCD and 
fragmentation is « 0.2 [5]. Second, protons and As do not seem to suffer from jet 
quenching. The measured nuclear modification factors Raa are close to 1 
very unlike those of mesons which show strong jet quenching Raa < 1 jHj- 



1 . 2. Hadronization of Bulk Matter 

In central heavy ion collisions we create a hot and dense fireball consisting of deconfmcd 
quarks and gluons. It is apparent that a hadronization picture that assumes a single 
parton in the vacuum has to fail in such an environment. Hadronization involving two 
or more partons can be described by higher twist fragmentation functions. However, 
nothing is known about these non-perturbative quantities. 

Instead we can look at the extreme case of a phase space densely populated with 
partons. In the vacuum, a single parton has to radiate gluons which subsequently 
split into quark-antiquark pairs. Hadronization can only start once sufficient quark- 
antiquark pairs are available to turn into valence quarks of hadrons. In a dense 
medium, enough partons are already there and could just recombine or coalesce into 
hadrons: three quarks into a baryon, a quark-antiquark pair into a meson. 

Recombination as a model of hadronization has a long history, see for some 
references. Recombination of beam fragments jH] has been experimentally confirmed 
by the observation of the leading particle effect in fixed target experiments I10| . The 
application to heavy ion collisions has been revived recently [111 El E3j and is a 
center of great interest since then. 

The yield of hadrons from a given parton system can be calculated starting from 
a convolution of Wigner functions (SJ- For a meson with valence (anti)quarks a and b 
we have 

d 3 N M f d3 R f d3 <ld 3 r Tlr / r P „ r P \ , . 

Here W a b is the 2-particle Wigner function for partons a, b and <I>m is the Wigner 
function of the meson. The sum runs over all possible parton quantum numbers. We 
note that recombination is dominated by the lowest state in the Fock expansion of 
the hadron, i.e. the valence quarks. Sea quarks and gluons do not seem to play a 
role, as can be inferred from the quark counting rule for elliptic flow (see below). For 
practical purposes the parton Wigner function is usually approximated by a product 
of single particle phase space distributions W a b = w a wi,. Several implementations of 
recombination have been discussed in the literature [5] El EH j see £Q for a review. 

In order to obtain an estimate when recombination is important as a hadronization 
mechanism, one can compare the yield of fragmentation and recombination starting 
from different parton spectra. The result of this competition is rather surprising. 
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Figure 1. Spectra of n°, p, Kq and A + A as a function of Pt at midrapidity 
in central Au+Au collisions at y/s = 200 GeV 0. Dashed lines are hadrons from 
recombination of the thermal phase, dotted line is pQCD with energy loss, solid 
line is the sum of both contributions. Data are from PHENIX (tt°, p) 1151 1^1 and 
STAR (K°, A + A) 0. 



Thermal, or at least exponential, parton spectra w ~ er p l T play a special 
role. Recombination on such a parton distribution leads to an exponential hadron 
distribution with the same slope since 

w a w b ~ e -*r/T e -(i-*)P/T = b -p/t_ (2) 

where x gives the momentum fraction of parton a. Therefore recombination is more 
effective than fragmentation on any thermalized parton ensemble. On the other hand 
one can show that a power-law parton spectrum favors fragmentation at least for large 
Pt, in accordance with perturbative QCD. 

Eq. Q does preserve momentum, but not energy. Therefore it can only be safely 
applied in a kinematic region where mass effects are small, i.e. for Pt M. On the 
other hand we know that pQCD fragmentation starts to dominate at very high Pt- 
The large jet quenching at RHIC suppresses the contribution from fragmentation, so 
that recombination effects can be observed at intermediate Pt- 

1.3. Recombination and RHIC data 

Calculations for Au+Au collisions at RHIC energies assume a thermalized system of 
constituent quarks with a temperature T around the phase transition temperature 
of the quark gluon plasma and strong radial flow. Recombination is applied to this 
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Pr(GeV) 

Figure 2. Schematic action of recombination and fragmentation on a transverse 
momentum spectrum of partons (solid line). To produce a meson with = 
6 GeV/c, the fragmentation process starts with a parton having much more 
transverse momentum. On the other hand, recombination can work with two 
coalescing partons having roughly 3 GeV/c each. 



phase. To describe the high- Pt spectrum of hadrons this has to be supplemented by 
a pQCD calculation using fragmentation and taking into account energy loss ^2 |3] . 
Alternatively one can use another model to describe the power-law tail of the Pt 
spectrum, e.g. recombination of hard partons or soft-hard recombination |12M14| . 

Fig. shows the Pt spectrum of 7r°, p, Kq and A + A in central Au+Au 
collisions obtained in 5] with a model including recombination and fragmentation. 
The agreement with available data is excellent for Pt > 2 GeV/c. We note that the 
hadron spectra exhibit an exponential shape up to about 4 GeV/c for mesons and 
up to about 6 GeV/c for baryons, where recombination of thermal quarks dominates. 
Above, the spectra follow a power-law and production is dominated by fragmentation. 

Recombination naturally leads to a ratio p/?r w 1. In the limit Pt —* oo 
where masses are negligible and kinematics is collinear, recombination of a thermal 
ensemble of partons exactly equals the statistical model of hadrons produced in 
thermal equilibrium. Therefore, recombination solves the RHIC puzzle of hadron 
enhancement by introducing a mechanism that can produce hadrons close to thermal 
equilibrium. It also pushes up the Pt of hadrons, see Fig. [21 The domination of 
recombination up to 6 GeV/c is also the reason why jet quenching seems to be absent 
in nuclear modification factors measured for baryons. Recombination predicts a sharp 
decrease in p/ir beyond 4 GeV/c which has not yet been seen in the data. 

l.^. Elliptic Flow 

Let us assume the parton phase exhibits elliptic flow v^ipr)- Recombination makes a 
prediction for elliptic flow of any hadron species after recombination [5] : 

v 2 (P T )=nvV(P T /n). (3) 

Here n is the number of valence quarks for the hadron. Hence recombination predicts 
another striking difference between mesons (scaling with n = 2) and baryons (scaling 
with n — 3). Different mesons should follow the same scaling law, even if the masses 
of the mesons are very different. 
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Figure 3. Elliptic flow 112 for 7r+, p, Kq and A as a function of Pt scaled 
by the number of valence quarks n vs Pr/n. The data follows a universal 
curve, impressively confirming the quark scaling law predicted be recombination. 
Deviations for the pions are discussed in the text. Data are taken from PHENIX 
(tt+, p) EH] and STAR (K°, A) Q. 



Fig- El shows the measured elliptic flow t>2 for several hadron species in a plot with 
scaled axes v^jn vs Pr/n. All data points (with exception of the pions) fall on one 
universal curve. This is a an impressive confirmation of the quark scaling rule and 
the entire recombination picture. The pions are shifted to lower Pt, because most 
pions in the detectors, even at intermediate Pt, are not from hadronization, but from 
secondary decays of hadrons. It was shown recently that inclusion of p recombination 
with subsequent decay of the p resonance into 2 pions can largely account for the 
shift ^7]. We note that the scaled elliptic flow in Figgis equivalent to the parton 
elliptic flow. This would be the first direct observation of a non-trivial observable in 
the parton phase. An an immediate consequence we conclude that strange quarks 
have the same elliptic flow as light quarks. 

The quark scaling law for elliptic flow is reminiscent of classical quark counting 
rules, like the famous 2/3 ratio of the 7r + N and N + N cross sections. Some decades 
ago, these counting rules were very important to convince us that quarks are indeed 
real degrees of freedom in high energy scatterings. Now we have discovered another 
quark counting rule at RHIC. We can rewrite JBJ as 

^ 2 M (2p) = 2 

vf(3p) 3 ( > 

where the superscripts M and B stand for the elliptic flow of a meson and a baryon 
respectively. This quark counting rule is a clear hint that quarks are relevant degrees 
of freedom at RHIC. Furthermore, elliptic flow is an observable to describe collective 
behavior. This implies that we created bulk matter with subhadronic degrees of 
freedom. Therefore, this discovery could be a very important step to prove the creation 
of a quark gluon plasma in Au+Au collisions at RHIC. 
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1.5. Further Remarks 

Recombination has been proven to be a valid description of hadron production in 
high energy heavy ion collisions. Hadron spectra, hadron ratios, nuclear modification 
factors and elliptic flow at intermediate momenta can be explained. The predictive 
power of recombination lies in the fact that a consistent description of all hadron 
species has been achieved with only one universal parametrization of the parton 
phase in terms of temperature, radial and elliptic flow and volume. Let us emphasize 
that the parton phase is input in the calculation, not a prediction. However, when 
hadron observables measured at RHIC are consistently described by one universal 
parametrization of the parton phase, confidence grows that the right parameters have 
been found. 

As mentioned above, recombination does not work at low Pt- This is not a 
failure of recombination itself, but a flaw of our simple description, that does not 
ensure energy conservation. Another problem that one encounters at low Pt is the 
requirement that entropy must not decrease during hadronization. Moreover, we know 
that hadrons at low Pt do not escape directly from hadronization, but they suffer 
multiple rescattering in the hadronic phase. Therefore, direct comparison with data 
at low Pt is difficult. Implementations of resonance production and decay into a 
recombination calculation can improve the description at low Pt ^3 E| • O n the 
other hand, measurements of resonances can be used to determine the strength of 
interactions in the hadronic phase 

1.6. Constituent Quarks? 

Let us take a closer look at the question what is actually recombining. From the 
quark counting rule for elliptic flow we can infer that only the valence structure of a 
hadron is relevant. Recombination involves only the lowest Fock state with sea quarks 
and gluons not being relevant. This seems to be at odds with our usual picture of a 
hadron. But this is not true. 

In high energy physics our picture of a hadron is mostly defined by parton 
distributions (PDFs) as seen in deep inelastic scattering. PDFs tell us that there 
is an infinite number of sea quarks and gluons in a hadron at a given Q 2 , where Q is 
a perturbative scale. Parton distributions refer to average parton configurations in a 
free hadron. 

We know that the momentum carried by sea quarks and gluons decreases if Q 
decreases. The evolution can not be continued quantitatively into the non-perturbative 
domain Q < 1 GeV/c. However, it is clear that for Q w Aqqd perturbative quarks and 
gluons are not longer appropriate degrees of freedom. Instead, one has to change to a 
constituent quark picture. In this picture quarks are dressed and acquire a dynamical 
mass and higher Fock states are suppressed due to the large quark mass. 

The constituent quark picture is not in contradiction with perturbative QCD 
because it lives in a domain where perturbation theory does not work. By applying 
non-perturbative methods in QCD it has been shown that there is a smooth transition 
between both pictures. Lattice QCD j20], instanton models [2] and solutions of the 
Dyson-Schwinger equation 22 provide consistent results of dynamical quark masses as 
a function of Q. These functions smoothly interpolate between current quark masses 
relevant for perturbative QCD at large Q and constituent quark masses at low Q. 
This running of the mass is also related to chiral symmetry breaking 22 . 
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At the phase transition with a temperature around 175 MeV, it can be ruled 
out that the deconfined phase is a perturbative quark gluon plasma. Instead, it is 
natural to have effective degrees of freedom with dynamical masses. A description of 
recombination that is linked to the chiral phase transition is not yet available, but 
should be a priority of future research. 

2. Hadron Correlations 

Recombination is thought to be the dominant mechanism of hadronization at 
intermediate Pt between 2 and 5 GeV/c. It has been argued that jet-like correlations 
of hadrons — observed in this momentum range — are an experimental fact that 
contradicts this view. This important challenge has only been addressed recently 2 . 

2.1. Associated Yields 

The STAR and PHENIX collaborations at RHIC have published data on dihadron 
correlations |24M25| . They detect a trigger hadron A in a certain Pt range and measure 
the yield of secondary hadrons B in another Pt range as a function of the relative 
angle A<j). The 2-hadron yield dNAB /d(Acf>) includes an uncorrelated background 
d(N 'aN b) I 'd(A<f>) . We define the background subtracted associated yield per trigger 
particle as 



It was found that the associated yield in Au+Au collisions at RHIC shows an 
enhancement above background around A(f> = which is approximately the same 
as that for p + p collisions. A smaller but broader enhancement can be seen in p + p 
around A(f> = ir, but it is absent in Au+Au collisions. 

In p + p collisions these signals can be interpreted as correlated emissions from 
jet cones. The near side peak (A(f> ss 0) is coming from hadrons formed in the same 
jet as the trigger hadron. The far side associated yield is due to hadrons from the jet 
recoiling against the jet of the trigger particle. It is washed out due to higher order 
corrections. The transverse momenta of both the trigger and associated particles are in 
the range of a few GeV/c. Therefore, there is no doubt that the correlations observed 
in p + p are due to jet fragmentation. 

On the other hand, near side correlations arc also observed in Au+Au collisions, 
while far side correlations seem to be severely dampened ESj- The absence of 
far side correlations can be understood as a consequence of the strong jet quenching 
in Au+Au at RHIC energies. However, the size of the near side correlations seem 
to be at odds with our understanding that a large recombination contribution in the 
relevant Pt range exists. 

2.2. Correlations from Recombination 

We have discussed in the first part how recombination was successfully applied to 
describe hadron spectra and elliptic flow starting from assumptions about the parton 
phase at hadronization. One crucial simplification that has been always implemented 




(5) 
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so far is a factorization of any n-parton Wigner function into a product of independent 
single parton distributions 

n 

H; „ \[n, (6) 

i=l 

By definition, this factorization does not permit any correlations between partons. 
Consequently, no hadron correlations can emerge via recombination. It has to be 
emphasized that the above factorization was chosen for simplicity and it was justified 
because single inclusive hadron spectra could be described very well. 

It has been proven in [2] that modifications of © including correlations 
between partons do lead to correlations between hadrons upon recombination without 
deteriorating the quality of the description of single hadron spectra. We will come 
back to this later. First, let us ask the legitimate question why there should be any 
correlations in the parton phase. The absence of correlations would indicate complete 
thermalization. However, the degree of equilibration at the time of recombination is 
not exactly known. 

Let us discuss one likely source of jet-like correlations in the medium. It has been 
found that the strong jet quenching observed at RHIC is due to energy loss of high-pr 
partons in the medium. The energy loss is estimated to be up to 14 GeV/fm for a 
10 GeV parton - This means that most jets apart from those close to the surface 
are completely stopped, dumping their energy and momentum into a cell of about 1 
fm 3 in the restframe of the medium. This results in a dramatic local heating, creating 
an ultra hot spot in the fireball. Moreover, the directional information of the jet is 
preserved. The subsequent evolution of the fireball might lead to partial diffusion of 
the hot spot. Nevertheless, residual correlations might still exist and mimic a jet, 
because the momentum of the jet has been transfered to the medium. If partons of 
such a hot spot are picked they would contribute to a jet-like correlation. Also a 
parton from a partially quenched jet together with a parton from the hot spot created 
can be correlated. The STAR Collaboration presented evidence that jet cones can not 
be seen as independent from the medium. Instead it is important that jets and the 
medium mutually influence each other |27|. 

2.3. Modeling Correlations 

For 2-meson production we have to consider the following possible scenarios: (1) Both 
mesons come from fragmentation, either of the same or of two different jets. In a 
straight forward notation we will call this process F-F. (2) Both mesons emerge from 
recombination of soft partons, denoted by SS-SS (S stands for soft). (3) While one 
meson recombines, the other emerges from the coalescence of a soft parton and a 
hard parton (from a jet). The latter is also called soft-hard coalescence, so we call 
this process SH-SS. (4) One meson comes from recombination while the other is from 
fragmentation (F-SS). (5) Both mesons can come from soft-hard recombination (SH- 
SH). (6) One meson fragments while the other is from soft-hard recombination (F-SH). 

All processes can produce both correlated and uncorrelated pairs. Here we are 
most interested in correlations coming from pure recombination (SS-SS). However, 
conventional jet correlations (coming from F-F) have to be taken into account when 
comparing to data. As an example for soft-hard correlations we will also estimate 
correlations from F-SH. These three processes are schematically depicted in Fig. 
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Figure 4. A schematic picture of correlated production of two mesons A and B 
through the processes F-F, F-SH and SS-SS. Thick lines denote soft partons. 



2.4- Dimeson Production from Soft Partons (SS-SS) 

We start with a simple extension of the correlation- free factorization ©. We want 

to assume that correlations are a small effect and restrict ourselves to 2-particlc 
correlations C,j , so that a 4-parton Wigner function can be written as 

W1234 « w 1 w 2 w 3 W4(l + ^Cij). (7) 

i<j 

The correlation functions C'ij between parton i and parton j are in principle of 
arbitrary shape, but we want to assume that they vary slowly with momentum and 
that they are only non- vanishing in a subvolume V c of the fireball. 

A picture of hot spots created by quenched jets motivates the following Gaussian 
ansatz in rapidity and azimuthal angle for the correlation functions: 

Cij = c So / oe -^-^) 3 /(20S) e -to-v,)'/&vS) (8) 

+ Cn S„ / 7re -(^»-^+^) 2 /(2^) g-to-w) 2 7(2i£) _ 

Here </>o j7r and yo,ir are the widths of the Gaussians in azimuth and rapidity, 
respectively. The two terms of the sum correspond to correlations initiated by an 
energetic parton (0 = 0) and its recoil partner (<fi = tt). Cq and c n give the strength 
of the near side and far side correlations, while the functions fo^{pTi,PTj) describe 
the transverse momentum dependence of the correlations. The functions So, ^ (<Xj, o~j) 
parametrize the spatial localization of the parton correlations on the hypersurface E. 
For simplicity we assume that So.w = 1, if 0$, Oj G V c and So,7r = otherwise. 

The 2-meson yield is given by a convolution of the partonic Wigner function W1234 
with the Wigner functions &a, $b of the mesons with an additional integration over 
the hadronization hypersurface S 

E a E b J d N ?JI =C ab [ da<f> A ®W 12 34 ®<pB- (9) 
d 6 P A d 6 P B J?, 

We restrict our discussion to near side correlations. To this end we neglect the 
second line in ©. We also work in a fixed Pt window and will therefore neglect 
the transverse momentum dependence by setting fo = 1, absorbing the normalization 
into the constant cq. With more data available in small Pt bins, the functional shape 
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of /o and /„. will be explored in the future. We assume Cq <C 1 which permits omitting 
quadratic terms like Cq or Coi>2- 

Using Boltzmann distributions for the single parton distributions Wi with 
temperature T, radial flow rapidity rjr, and a boost-invariant hypersurface S at proper 
time r, Eqs. ©, J7J and JSJ lead to 

P wp m dP m d =( 1 + 2d '° + Uoe- 1 **?'™') (10) 

Pta dPrA d(f>A dy a Ptbo-Ptb d<PB dyB \ ' 

xh A {P TA ) (1 + 2v 2A (P TA ) cos(20 A )) h B {P TB ) (1 + 2v 2B {P TB ) cos(20 B )) . 

V2i is the elliptic flow coefficient for hadron i and A(f> — \<j) A — <j) B \ is the relative 
azimuthal angle. The factor Q = 4 in front of the Gaussian term is called the 
amplification factor. We have introduced the short notation 



where rrij are the masses of the recombining (anti)quarks, Mj is the transverse mass 
of the meson, At is the transverse area of the fireball and Ci is the degeneracy factor 
for meson i. The integration over the correlation volume has been absorbed into the 
correlation strength leading to the new parameter cq « cqV c /(tAt). 
The single meson spectra take the simple form 

d 3 N 

p -775 ZT^ = ^( p n)(l + co) (1 + 2« 2i (P ri ) cos(2^)) . (12) 
P T idP T id4>idyi 

We note that the only modification to the single meson spectrum — compared to the 
case without correlations [H] — consists of a moderate rescaling 1 — > 1 + cq. For small 
numerical values of Co this can be easily absorbed in the overall normalization, so that 
a good description of single particle spectra in the recombination formalism can still 
be achieved with a consistent set of parameters for the parton phase. 



2.5. Amplification of Correlations 

The factor Q = 4 implies an enhancement of the correlations on the hadron side 
compared to the parton side. The effect is essentially the same as in the amplification 
of elliptic flow by the number n of valence quarks in the hadron. In the case of 2- 
parton correlations, Q counts the number of possible correlated pairs between the n A 
(anti)quarks of meson A and the n B (anti)quarks of meson B. Since we work in a weak 
correlation limit where quadratic terms are suppressed, we do not take into account 
more than one pair of correlated partons. Apparently 

Q = n A n B , (13) 

thus Q = 6 for a meson-baryon pair and Q = 9 for a baryon-baryon pair. The term 
2co in Eq. (|10[) comes from the possible correlation of the quark-antiquark pair inside 
of one of the two mesons. 

We can now study the associated yield defined in Eq. |JSj, in a certain kinematic 
window, by subtracting the uncorrelated yield N A N B , determined by (|12f> . from the 
total dimeson yield given in (|10|) . The terms proportional to 1 + 26o cancel and hence 

2ttN a Y ab (A<Z>) = Qc a e~ iA ^ 2 /^ 2 N A N B . (14) 
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The Ni are single particle yields in the kinematic window of the trigger meson or 
associated meson 



Ni = 2?r / dy 



: J dPrPrhi(Pri). (15) 

Correlations in rapidity have not been studied here. Since we integrate over the 
rapidities any residual effects can be absorbed into the constant cq. We note that the 
result in l|14|) also holds for meson-baryon and baryon-baryon correlations with the 
appropriate amplification factor Q. 

It is also interesting to know the size of the subtracted background. It is given 

by 

WA ^^f = NaNb (1 + 250 + 2 ^ 2B cos(2A0)) . (16) 

V2i is the elliptic flow of meson i averaged over the respective Pr window. We note 
that the background already has a non-trivial dependence on the relative azimuthal 
angle coming from the the elliptic flow of the mesons. 

2. 6. Correlations from Fragmentation 

Before comparing to data we have to find a way to describe the presence of 
correlations from the F-F process. This requires a study in the framework of dihadron 
fragmentation functions. These functions are under investigation, but their functional 
form is still poorly known [2H|- Therefore, we choose a simple model based on single 
hadron fragmentation functions. After a parton a with momentum p has fragmented 
into a hadron A with momentum Pa = zap, described by a fragmentation function 
Do,/a(.za)i we assume that the production of the second hadron B with momentum Pb 
is given by a fragmentation function D cl / B (zb) where Pb = zb(1 — za)p- This model 
has clear limitations, but it should give a rough estimate of the transverse momentum 
dependencies. 

Furthermore we introduce a Gaussian distribution of the relative azimuthal angle 
and rapidity of the two hadrons A and B. It can be shown that 

d{N A Y A B{&4>)) o_r/o_^-l/2„-(A<*) 2 /(2^) P dz A 



27r/(27r^)-i/2 e -CA0) 2 /(20o 2 ) £ I 

„ J Zn 



dPr A dPTB a-Jzo z a{1-z a ) 

9aip) — E a d 3 N a /dp is the spectrum of the fragmenting parton a and I is the result 
of integrating the rapidities of the two hadrons in the kinematic window. O is the 
width of the azimuthal correlation and the limits for za are given by Zq = 2PtaI\[s 
and z\ = Pta/(Pta + Ptb)- We also included a possible energy loss AE of parton a 
before fragmentation. 

The fragmentation from a single parton provides correlated hadrons. The 
uncorrelated background can be easily estimated from fragmentation of two hadrons 
from two independent jets. 

2. 7. Correlations from Soft-Hard Processes 

Among the possible 2-hadron production mechanisms involving soft-hard recombina- 
tion, we choose the F-SH process as an example for our discussion. Soft-hard recom- 
bination is not part of the original recombination model in [111 [S] , but it is taken into 
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account by other groups |12M14) . The relative importance of soft-hard recombination 
is not clear, since both approaches fit the experimental data well. It was argued that 
soft-hard processes are a good source of hadron correlations. We want to shed some 
light on this discussion. 

We use the following model for dimeson production in the F-SH channel. Suppose 
a hard parton a with momentum p a produces an additional bb pair through branching, 
where b is any light quark flavor. The pair ab hadronizes into a meson, while b can 
pick up another (anti)quark c from the surrounding medium to coalesce into a meson 
B. The production of A can be described by a fragmentation function D a ^(zA) while 
we want to assume here that parton b remains with momentum (1 — ZA)p a - 

With a Boltzmann distribution for parton c it follows that 

d(N A Y AB (A(j))) ...SCbMtB _ , n/o_.,2n-1/2.,-(A«A) 2 /(2^) r ( ^TB SUlh T]T " 



dP T AdPTB 



xg ,i ^±gg^ (18) 



with p a = Pta + Ptb/2 + AE. v < 1 is a normalization constant arising from the 
restriction of the S integration to a subspace given by the intersection with a jet cone. 



2.8. Numerical Results 

For the numerical evaluation we use the set of parameters found in [5] to fit the single 
hadron spectra and elliptic flow measured at RHIC. We the minijet calculation in 
and KKP fragmentation functions j3U]. We choose the windows 2.5 GeV/c< Pta <4.0 
GeV/c for trigger particles, 1.7 GeV/c< Ptb <2.5 GeV/c for associated particles and 
\ua\, \ub\ < 0.35 as in the recent analysis of the PHENIX experiment We fix 

the azimuthal correlation width to be 4>o — 0.2 in all channels, in rough agreement 
with the experimental data. In order to have a measure of the absolute strength of 
the correlation we follow (23 and integrate Yab over the near-side peak 

/■0.94 

YXT= / d(A$)F AB (A$). (19) 
Jo 

Fig. El shows the integrated associated yield of hadrons for the case that the 
trigger is a baryon (proton or antiproton) and a meson (pion or kaon) for different 
centralities. We consider the following cases: (i) F-F, i.e. fragmentation only. Our 
results are in rough agreement with the data without any adjustable parameter in 
our model. This means that our model for double fragmentation would be suitable to 
describe correlations in p + p. (ii) F-F and SS-SS with cq — 0. Keeping fragmentation 
and turning on recombination without correlations dilutes the signal as expected. The 
effect is particularly strong for baryon triggers, (iii) F-F and SS-SS with co = 0.08, 
and (iv) F-F and SS-SS with co = 0.08 x 100/iV pa rt- We remember that Co contains 
a space-time integral which can depend on centrality. Therefore we test two different 
scenarios, cq constant and Co scaling like ~ N~^ rt . These cases correspond to the 
correlation volume increasing, V c ~ -/V part , and V c being held constant respectively. 
The first case leads to a strong rise of Y A *g e with increasing centrality that is not 
observed in the data. The scenario of approximately constant correlation volume 
V c leads to a nearly constant Y2°g c which is in qualitative agreement with PHENIX 
measurements |25|. (v) We also consider the case F-SH (for pions only). A choice of 
it = 0.5 for the parameter leads to sizable correlations. However, the yield of hadrons 
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Figure 5. The associated yield Yf£ e for baryon triggers (right panels) and 
meson triggers (left panels) as a function of iVp ar t. Squares are F-F only (i), 
diamonds are F-F and SS-SS with co = (ii), triangles are F-F and SS-SS with 
c = 0.08 (iii), circles are F-FF and SS-SS with c = 0.08 X lOO/Npart (iv), stars 
are F-SH with v = 0.5 [pions only](v). Insert: associated yield Y^s a s a function 
of A0 before background subtraction from F-F and SS-SS at impact parameter 
b = 8 fm. 



from the F-SH channel is small compared to that from F-F and SS-SS. This is due to 
the small contribution from soft-hard recombination in the Duke parametrization as 
already discussed in [5]. 

We note that a realistic model of correlations from fragmentation and 
recombination can reproduce associated yields in qualitative agreement with the data. 
The parameters of the parton phase can be chosen consistently so that the excellent 
fits of single hadron spectra and elliptic flow are preserved. Correlations on the parton 
level need only to be 8% in order to reproduce the data. Soft-hard recombination 
naturally leads to correlations between hadrons, but this hadronization channel is not 
necessary to obtain a good description of the data. More details can be found in 



Hadron Correlations from Recombination and Fragmentation 



14 



3. Conclusions 

I have presented evidence that hadronization in heavy ion collisions is dominated 
by recombination from a thermalized parton phase. A theoretical description of 
recombination for Pt > 2 GeV/c is available and naturally explains the enhanced 
baryon production at RHIC. The quark counting rule for elliptic flow, impressively 
confirmed by experiment, might prove to be an important cornerstone to make the 
case for the quark gluon plasma. 

I also emphasized that our description of the dynamics of the recombination 
process is still insufficient, which hampers progress in our understanding of low Pt 
hadron production. In particular, we do not have a quantitative theory that connects 
recombination to the chiral phase transition. 

On the other hand, it has been shown that experimental findings of jet-like 
hadron correlations at intermediate Pt are not in contradiction with the existence 
of recombination. Correlations among partons can be introduced into the existing 
formalism and naturally lead to correlations among hadrons. Jet-like correlations in 
the medium can be created by jets that are quenched by the medium and leave a hot 
spot. 
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